Surface density of states of s±-wave Cooper pairs in a two-band model 
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We calculate the surface density of state (SDOS) of s±-wave Cooper pair in two-band super- 
conductor model, where gap functions have different signs between two bands. We find that the 
Andreev bound state appears at surface due to the sign change in the gap function in the interband 
quasiparticle scattering. However, we do not obtain the zero-energy peak of SDOS in contrast to 
the d-wave case. The tunneling spectroscopy of s±-wave is much more complex as compared to the 
d-wave case realized in high-T c cuprates. 
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I. INTRODUCTION 

Recent discovery of superconductivity in the iron based 
LaFeAsOi-zFa; with T c = 26K— has aroused great inter- 
ests as a class of non-cuprate compound. In the iron- 
based family, various compounds exhibit superconduc- 
tivity with T c now exceeding 55 K. Superconductivity 
has also been found in iron-based materials with dif- 
ferent layered structures that include BaFe2As2^ and 
FeSe 3 . Local spin-density calculations for LaFeAsO have 
shown that the system is around the border between mag- 
netic and nonmagnetic states, with a tendency toward 
antifcrromagnetism4i£ It has also been pointed out that 
the electron-phonon coupling in this material is too weak 
to account for T c = 26K. 6,7 Based on the first principles 
calculation, minimum five-band model to describe the 
iron-based superconductor has been proposed 8 -. Using 
this five-band model, pairing symmetry has been calcu- 
lated based on the random phase approximation (RPA) 8 . 
The resulting gap function does not have nodes on the 
Fermi surface while it has a sign change between Fermi 
surfaces. Now, it is called an s±-wave pairing 9 "^. There 
have also been relevant theoretical predictions which sup- 
port the realization of the s±-wave mode l 11 ' 12 ' 13 . 

In order to elucidate the energy-gap structure of these 
s±-wave superconductors, experiments based on stan- 
dard technique, e.g., NM R 14 ' 15 , specific heat£, penetra- 
tion dept h 16 ' 17 and quasiparticle tunneling spectroscopy 
have starte d 18 ' 19 ' 20 ' 21 . It is a very challenging issue to 
clarify the superconducting profile of s±-wave supercon- 
ductors. Since the internal phase degree of freedom ex- 
ists in the gap function of s±-wave pairing, it is natural 
to expect phase sensitive phenomena realized in high- 
T c cuprate a 22 ' 23 ! 24 ' 25 . As shown in the study of high-T c 
cuprates, the mid gap Andreev bound state (MABS) is 
formed at the surface due to the internal phase effect, 
where a quasiparticle feels a different sign of the gap 
function depending on the direction of their motions. 
The presence of the MABS produces zero-energy peak 
(ZEP) of the surface density of states (SDOS) and has 
been observed as a zero-bias conductance peak (ZBCP) 
in tunneling spectroscopy up to no w 26 ' 27 ' 28 ' 29 . It is an ur- 



gent topic to reveal whether MABS exists in the s±-wave 
pairing or not. 

Several theories of surface or interface profiles 
about s-t-wave pairing have been presented very 
recently^ 3 !^^^^^^ 7 -. However, SDOS of two-band 
superconductors with s± model has not been understood 
yet. The presence or absence of MABS have not been 
resolved yet. Furthermore, character of the nonzero in- 
ner gap Andreev bound state (ABS) has not been clar- 
ified. To reply to these issues, in the present paper, we 
employ a simple two-band tight-binding model with s±- 
wave as a prototype of iron-based pnictides, and calcu- 
late SDOS for the [100] and [110] oriented interfaces us- 
ing the t-matrix method 3 ^. A merit of our calculation 
is that details of the band structure and band mixing 
can be microscopically taken into account. We find that 
ABS with nonzero energy is formed at the surface due 
to the interband quasiparticle scattering, through which 
gap functions change sign. However, there is no ZEP in 
SDOS in contrast to the case of cZ-wave pairing realized 
in high-r c cuprate o 26 ' 27 ' 2 ^. 

II. MODEL AND FORMULATION 

We start with a two-band tight-binding model on a 
square lattice. Since there has not been an explicit study 
about the SDOS of the s± model, we choose the energy 
dispersion of two orbitals simply supposing d xz and d yz 
orbitals or p x and p y orbitals. Hereafter, we define index 
1 in matrix form as d xz {p x ) orbital and index 2 as d yz 
(p y ) orbital. X and Y axes are rotated by 45 degrees 
from x-y, where x and y denote the axes in a unit cell as 
shown in Fig. [T] 

First, we discuss the normal state. Tight-binding 
Hamiltonian is given in the form 

#o = tvJ^LaCjw C 1 ) 

ij iiv u 

where fj M '„ is a hopping integral from the v-Xh orbital 
on the j-th site to the /i-th orbital on the i-th site, cj 
creates an electron with spin a on the ^-th orbital at 
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FIG. 1: (Color online) Two-dimensional square lattice with 
nearest-neighbor hopping t (thin solid line) and next-nearest- 
neighbor hopping t' (thin dotted line). The [100] and [110] 
oriented surfaces are constructed by inserting four infinite po- 
tential barriers illustrated with thick solid lines and thick dot- 
ted lines, respectively 



site i. As shown in Fig. [TJ we take the nearest-neighbor 
hopping t and the next-nearest-neighbor hopping t'. The 
band filling n is defined as the number of electrons per 
number of sites (e.g., n — 1 for half filling). The Hamil- 
tonian in the Fourier transformed representation is given 
as 
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where the 2x2 matrix i°(k) is denoted by 



f°(fe) 



— t cos k x 
2t' sin k r . sin k„ 



sin k x sin ky 
—t cos k v 



(2) 



(3) 



Hereafter, we take t and the lattice constant a as the 
units for energy and length, respectively, i (k) can be 
diagonalized to £„(&), which corresponds to the energy 
of band a 



(4) 



where U(k) is a 2 x 2 unitary matrix. Fermi surfaces 
consist of two parts near the half filling as shown in Fig. 
[21 We define the band which forms inner (outer) Fermi 
surface as band — (+). 

In a two-band model, gap function generally forms a 
2x2 matrix. The gap function A^„(fe) in the orbital rep- 
resentation is transformed to the gap function A a t(k) in 
the band representation using the unitary matrix U(k), 



A ab (k) = ]T u; a (k)u; b (-k)A^(k). 



(5) 



Here, we neglect the frequency u> dependence of the gap 
function and assume that the gap function in the band 
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FIG. 2: Outer (inner) Fermi surfaces with gap function A + (_\ 
consist of band +(— ). 



representation is diagonal, A a (fe) = A aa (fe). Gap func- 
tion of band — (inner Fermi surface) and band + (outer 
Fermi surface) are denoted by A_ and A + , respectively 
as shown in Fig. [5] Then, the gap function in the orbital 
representation is obtained as 



A^fc) 



U^{k)U va {-k)A a {k). 



(6) 



In the case that the hopping integral is given by a 
real number, the relation e°„(fc) = £°^ v (—k) is satisfied. 
Then, we take following relation 



£W*0 = u*(-k). 



(7) 



Using the above gap function, bulk Green's function 
G(w, k) in the superconducting state is given by a 4 x 4 
Nambu representation as follows 



G(w,fe) 



e°(k) -pi 
At(fe) 



A(fe) 
-i°(k)+fx 



(8) 



with chemical potential [i. In the actual numerical 
calculation, we replace u> by lj + ij with small real 
number 7 to avoid divergence of the integral. Lo- 
cal density of states (LDOS) in the bulk is obtained 
by — l/(Nn) J^ fe ;=1 2 ^Gu(ui, k), where N denotes k- 
point meshes. The Green's function of the inhomoge- 
neous system including surface G s (u>,r,r') is calculated 
by G(u>,r), which is the Fourier transformed form of 
G(cu,k). As shown in Fig. [TJ we insert the infinite po- 
tential barrier Z(r) in four-atomic layers parallel to the 
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FIG. 3: (Color online) SDOS of the [110] and [100] ori- 
ented surfaces and LDOS in bulk are depicted with solid line, 
dashed-dotted line and dotted line, respectively for n — 1, 



t' = 0.1, and A- 
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These values are normalized by the values in the normal state. 

surface in the actual calculation. G s (ui, r, r') is given by 

G s (uj,r,r') = G(uj,r-r') 

dr"G(uo, r - r")Z(r")T 3 G%Lo, r", i$) 



where f denotes the Pauli matrix in charge space. 
We note that G s breaks translational symme- 
try. Using the G s (uj,r,r'), we obtain SDOS by 
— l/7r^ z=1 2 ImGf ; (w, r s , r s ), where r s denotes the 
location of the surface. Throughout this study, we take 
N = 4096 x 4096 fc-point meshes and 7 = 0.003. 



III. RESULT 
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FIG. 4: (Color online) Plots of LDOS similar to Fig. [3] for 
n = 1, t' = 0.1, and A_ = A+ = 0.1 (s-wave). 
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FIG. 5: (Color online) Plots of LDOS similar to Fig. [3] for 
n = 1, t' = 0.1, and A- = — A+ = —0.1 (s±-wave). 



In the following, we focus on LDOS at the surface, 
i.e., SDOS and bulk. First we show the result of the 
d 2; 2_j / 2-wave case for n = 1, t' = 0.1 in Fig. [3] where 
the gap function of band +(— ) is chosen as A + (_) = 
0.05(cos h x cos /by). Throughout the present study, 
LDOS is normalized to that of the value in the normal 
state at w = 0. We see a sharp ZEP of SDOS in [110] 
oriented surface and V-shaped LDOS in the bulk, which 
are consistent with the case of the single band d x 2_ y 2- 
wav o 26 i 27 ' 28 . The origin of the sharp ZEP is the sign 
change in the gap function felt by quasiparticles scattered 
at the surface, where the momentum of the quasiparticles 
parallel to the [110] surface is conserve d 26 i 27 ' 28 . 

As a reference, we show the result of s-wave pairing 
in Fig. GQ where n = 1, t' = 0.1 and A_ = A + = 0.1. 
The line shapes of SDOS of the [100] and [110] oriented 
surfaces, and LDOS of bulk are almost identical. This 
behavior is robust irrespectively of the band structure as 
far as the relation A_ = A + is satisfied. 



Next, we move to the s±-wave case. The corresponding 
results for s±-wave with n = 1, t' = 0.1 are shown in Fig. 
where we choose A + (_) = +(— )0.1. We see that two 
sharp peaks within the bulk energy gap appear in the 
[110] oriented SDOS^i. The clear difference from Fig.3 
is that there is no ZEP. For the [100] oriented surface, 
the value of the corresponding SDOS within the gap is 
almost constant with nonzero value. 

In order to clarify the origin of the two peaks in 
SDOS of the [110] oriented surface, we show the ky- 
resolved SDOS of the [110] surface in Fig. [6] where 
ky denotes the momentum parallel to the [110] sur- 
face. The fc-^-resolved SDOS is enhanced at ky ~ 
±\/27r/4, ±3\/27r/4, which correspond to (k x ,k y ) — 
(0, ±71-/2), (±71-/2, 0), (±tt, ±7t/2), (±7r/2, ±tt) at the orig- 
inal Fermi surface as shown in the inset of Fig. [5] At 
these points, angle resolved SDOS has a large value. Fur- 
thermore, quasiparticles feel a different sign of the gap 
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FIG. 6: Contour plot of the fcy-resolved LDOS(SDOS) in the 
[110] oriented surface for n = 1, t = 0.1, and A_ = — A+ = 
—0.1 (s±-wave), where ky = ±\/2/4 and ky = ±3\/2/4 are 
depicted as dotted lines. In the inset circles denote points 
on the Fermi surface with ky — ±V2/4 and ky = ±3v2/4, 
which mainly contribute to the LDOS. 



function through the scattering between inner and outer 
bands, which brings about the ABS. The large momen- 
tum change in quasiparticlcs is automatically induced by 
infinite potential barriers inserted at the surface due to 
normal (backward) reflection. 

On the other hand, the fc^-resolved SDOS in the [100] 
surface (not shown) is enhanced at k y — ±7r/2 within 
the gap. Since scattering of the quasiparticle preserving 
k y = ±7r/2 occurs between the inner and outer Fermi 
surfaces, ABS appears within the bulk energy gap, which 
makes residual LDOS, shown as the dashed-dotted line 
in Fig. El 

Although the sign change in the gap function does not 
produce the ZEP as in the case of unconventional super- 
conductors such as d- or p-wave pairing, the sign change 
in the gap function felt by quasiparticle enhances the 
magnitude of the inner gap LDOS. In a certain case, it 
produces sharp peaks, shown as a solid line in Fig. [5] In 
order to confirm whether the above behaviors are robust 
or not, we change the shape of the Fermi surface by con- 
trolling the value of t' for n = 1, A_ = — A + = —0.1. 
As shown in Fig. [7J the positions of the two peaks in 
SDOS for the [110] oriented surface (solid line) move to- 
ward that of bulk LDOS (dotted line) as the value of 
t' increases. Thus, we find that the positions of peaks 
of the [110] oriented SDOS are sensitive to the shape of 
the Fermi surface. The resulting peak positions are rele- 
vant to the relative position between the outer and inner 
Fermi surface. 

It is also interesting to clarify how the above results are 
influenced by changing the value of A_. In Fig. we 
focus on the A_ dependence of SDOS and bulk LDOS 
for A + = 0.1, n = 1 and t' = 0.1. Two-gap structure 
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FIG. 7: (Color online) Plots of LDOS similar to Fig. ES for 
n = I, A_ = -A+ = -0.1, t' = 0.3 (left panel) and t' = 0.5 
(right panel). 
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FIG. 8: (Color online) Plots of LDOS similar to Fig. [3]for n = 
1, t' = 0.1, A_ = -0.5A+ (top left panel), A_ = -0.1A+ 
(top right panel), A_ = 0.5A+ (bottom left panel), and A_ = 
0.1 A+ (bottom right panel). 



appears in bulk LDOS (dotted lines). It is very clear that 
the resulting bulk LDOS is insensitive to the sign of A-_ 
by comparing the A„ = 0.5A+ (A_ = 0.1 A + ) case with 
the A_ = -0.5A+ (A_ = -0.1A+) case. As far as we 
are looking at bulk LDOS, there is no difference between 
s±-wave and s-wave. The internal phase degree of the 
gap function does not appear in the bulk LDOS. On the 
other hand, sharp peaks of ABS appear only for SDOS 
of the [110] oriented surface with A_ = — 0.5A+ and 
A_ = — 0.1A + . The position of the sharp peaks moves 
toward to — with the decrease of the magnitude of A_. 
At the same time the height of the peaks is reduced (solid 
lines in the upper two panels). On the other hand, SDOS 
for the [100] oriented surface does not have clear peaks 
as compared to that for the [110] oriented surface. 

Finally, we show the result of n — 1.2 for t 1 = 0.1 and 
A_ = -A+ = -0.1 in Fig. [9] In this case, SDOS for 
the [110] oriented surface has two peaks at ui = ±A + in 
addition to the two inner gap peaks. This is understood 
by using the ky -resolved SDOS in Fig. [TU1 We see that 
ABS within gap vanishes for — 2v / 2~tt/3 < ky < —V2n/3 
and \/2V/3 < ky < 2^/2i:/2> since interband pair scat- 
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FIG. 9: (Color online) Plots of LDOS similar to Fig. [3] for 
n = 1.2, t' = 0.1, A_ = -A+ = -0.1. 
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FIG. 10: Contour plot of the fcy-resolved LDOS(SDOS) in 
the [110] oriented surface for n = 1.2, t' = 0.1, and A_ = 
-A+ = -0.1, where k Y = ±V2/3 and k Y = ±2^2/3 are 
depicted as dotted lines. In the inset we see that there is no 
interband pair scattering for — 2v2tt/3 < ky < — v2~7r/3 and 
V2n/3 <k Y < 2V2tt/3. 



are completely different from the d-wave gap function 
in high-T c cuprates where ZBCP appears for the [110] 
orientatio n 26 ! 27 ' 28 ! 38 . One of the big difference from the 
d-w&ve case is that quasiparticles do not feel the sign 
change in the gap function as far as it is scattered within 
the same band in s±-wave pairing. There are always 
both intraband and interband pair scattering. The pres- 
ence of intraband scattering without a sign change may 
prohibit the generation of MABS and ZEP of SDOS. In 
the present study, we employ a simple two-band model. 
We admit that five orbitals are needed to describe the 
superconductivity of iron-based superconductors. Our fi- 
nal goal is to establish a theory of tunneling spectroscopy 
/ surface density of state taking into account five bands. 
However, up to now, there has not been fully microscopic 
theory of surface density of states of multi-band super- 
conducting systems even in two-bands cases. To under- 
stand the essence of the interference effects originating 
from the existence of the multiband, it is reasonable to 
start with the two-band model. Thus, in the present pa- 
per, we have chosen the two-band model for the first step 
as a prototype of multiband model. In the near future, 
we will report the results based on the more realistic five- 
band model. 

Finally, we comment about the relevance of the present 
s±-wave model in two-band systems and two-band model 
in non-centrosymmetric superconductors^. Recently, 
there are several studies about surface density of states 
of non-centrosummetric superconductors. The presence 
of the Rashba spin-orbit coupling induces Fermi surface 
splitting, and a similar situation in the present two- 
band model seems to be realized. Due to the presence 
of Rashba interaction, the spatial inversion symmetry is 
broken in these systems. Then spin-singlet s-wave and 
spin-triplet p-wave pairing can mix each other. If the 
magnitude of the p-wave component is larger than that of 
s-wave one, ABS exists and MABS is possible for the per- 
pendicular injection of the quasiparticlc. The resulting 
ABS can be regarded as helical edge modes and carry spin 
current. The direction of the current flow corresponding 
to each Kramers doublet is opposite. On the other hand, 
the profile of the ABS in the present two-band s±-wave 
model is very different. In the present case, there is no 
spin current and spin degeneracy remains. 



tering is prohibited due to the absence of outer Fermi 
surface as shown in the inset of Fig. [TO] Thus, in this ky 
region, the angle-resolved LDOS is an independent sum- 
mation of LDOS in the outer and inner Fermi surface. 
Then, the resulting angular averaged LDOS has peaks at 
w = ±A+. 

Summarizing the above results, although the quasipar- 
ticle feels a sign change of the gap function for fixed 
ky {ky) in the reflection process at the surface, ZEP 
does not appear in the SDOS of s±-wave pairing. Thus, 
ZBCP in tunneling spectroscopy appears neither the 
[100] nor the [110] oriented junctions. These features 



IV. CONCLUSIONS 

In this paper, we have calculated the [100] and the [110] 
oriented SDOSs for s±-pairing in a two-band model by 
changing the shape of the Fermi surface and the band fill- 
ing. It has been revealed that the inner gap sharp peaks 
appear for SDOS in the [110] oriented surface. These 
peaks originate from the ABS caused by the interband 
scattering of quasiparticles, through which gap functions 
change sign. Such sharp peaks do not appear in the s- 
wave case, where there is no sign change in the gap func- 
tion between the two bands. It is also noted that the 
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resulting SDOS of s± model does not have ZEP. This 
means that the tunneling spectroscopy of s± supercon- 
ducting state is much more complex as compared to the 
d-wave case realized in high-T c cuprates. 

Up to now, there has been experimental reports about 
tunneling spectroscopy. The experimental line shapes 
of tunneling conductance are distributed including gap 
struct ure o 15 i 16 ' 17 i 21 and ZE P 18 ' 19 ' 20 . However, the exper- 
imental condition has not been clarified yet up to now. 
In the light of the study of high T c cuprat o 40 ' 41 , it has 
been revealed that well-oriented surface or well-oriented 
interface with low transparency junctions are needed to 
compare the surface density of states with the actual tun- 
neling conductance^ 2 -. We hope tunneling spectroscopy 
of well-oriented surface or well-controlled junctions with 
low transparency will be attainable in the present iron- 
based superconductors by the progress of microfabrica- 
tion technique. 

There are several interesting future problems. In the 
present paper, we have solved the Green's function in 
tight-binding model. It is possible to solve the Bogoli- 
ubov de-Gennes equation in the lattice model. The study 



along this direction is useful to elucidate interference ef- 
fect much more in detai l 43 ' 44 . Josephson effect in the 
s±-wave superconductor may be fascinating since we can 
detect internal phase effec t 45 ' 46 . It is interesting to clar- 
ify the possible existence of nonmonotonic temperature 
dependence in high-T c cuprate junctions 47,48 . Proxim- 
ity effect in s±-wave superconductors is also an inter- 
esting topic. Through the study of the proximity effect 
in unconventional superconductors 4 ^, the odd-frequency 
pairing amplitude has a crucial role to characterize the 
bound stated. It is a challenging issue to clarify the in- 
duced odd-frequency pairing near the present two-band 
model. 
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